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Abstract .  The  problem  of  intact  loading  cf  a  canti¬ 
lever  bean  which  exhibits  either  strain-hardening  or  strain- 
rate  sensitivity  is  considered .  Cn  the  basis  cf  an  assump¬ 
tion  concerning  the  node  cf  defcrnaticr.  of  the  bear,  tie 
the  bear  is  analysed.  The  influence  cf  strain- 
ana  strain-rate  sensitivity  on  the  remanent  de- 

In  addition,  numerical  res -alts  are 
presented  for  the  predicted  plastic  strains  and  strain-rates 
within  the  bean,  Finally,  directions  in  which  an  extension 
cf  the  analysis  is  desirable,  are  indicated. 


notion  of 
hardening 
formation  is  calculated. 


The  results  presented  in  this  paper  were  obtained  in 
the  course  of  research  sponsored  by  the  Iff ice  of 
F&val  Research  'under  Contract  Ucnr- 562(10)  with  Brown 
Research  Assistant,  Division  of  Applied  Mathematics, 
Brown  University,  Providence  12,  Rhode  Island. 
Professor  of  engineering ,  Brown  University. 
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distance  from  neutral  axis  to  extreme  fibre 
coefficient  in  c  -  s  relation 
coefficient  in  moment-curvature  relation 
Young's  modulus 

5  2 rj  +  p2  •*.  (1  --  (1  +  p)2  -  1  -  -In  (1 

A  a  +  i)2  -  i) 


w>c  di._ 


i  - 


mass  at  tip  of  beam 

=  -  1 

P 

moment  of  inertia  cf  cross-section  of  beam 

=  ^  =  ratio 

nvt  be  am  ma  s  s 

length  of  beam 

(1,  mass  per  unit  length  of  beam. 

(2)  as  subscript,  denotes  quantities  pertaining 
root  cf  the  beam 


bending  moment 

static  limit  moment  of  beam 

ratio  7- — r —  of  slopes  cf  elastic  ana  clastic  c-arts 
(do/  de  )p£ 

of  d-e  curve 

exponent  in  stress-strain-rate  relation 

shearing  force  in  beam 

ratio  _ initial  kinetic  energy 

maximum  elastic  strain  energy 

ratio  IV M0 

initial  value  of  f 
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6 0  static  yield  stress 

cp  slope  of  beam  =  ~ 

dx 

•  dot  denotes  differentiation  with  respect  to  time 

1.  Introduction.  We  consider  a  uniform  cantilever 
beam  with  attached  tip  mass  moving  initially  with  velocity  V 
normal  to  its  length.  At  time  t=G  the  root  of  the  beam  is 
instantaneously  brought  to  rest  by  contact  with  a  rigid  stop, 
(See  Fig.  1).  The  problem  is  to  determine  the  subsequent 
motion  of  the  beam  and,  particularly,  its  permanent  plastic 
deformation.  Interest  in  this  problem  stems  from  the  fact 
that  it  represents  a  case  in  which  the  plastic  deformation 
of  beams  under  impact  loading  can  be  conveniently  investi¬ 
gated  experimentally. 

Lee  and  Symonds  [l]  showed  that  a  comparatively  simple 
theoretical  analysis  of  problems  of  dynamic  loading  of  beams 
is  possible  if  it  is  assumed  that  the  beam  is  made  of  a  rigid- 
ideally-plastic  material.  This  rigid-ideally- plastic  analysis 
was  further  developed  by  Symonds  and  Leth  [2]  and  by  Men tel 
[3] •  In  addition,  experimental  investigations  of  the  pro¬ 
blem  have  been  carried  out  at  Brown  University,  and  the  re¬ 
sults  have  been  reported  by  Symonds,  Green  and  hen  tel  [4], 
and  by  Mentel  [3].  Large  quantitative  discrepancies  were 
found  between  theory  and  experiment.  Mentel  [3]  attributed 
the  discrepancy  to  strain-hardening  and  strain-rate  effects 
which  are  neglected  in  the  rigid-ideally  plastic  analysis. 
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Impact  tests  made  by  Parkes  [ 5l  showed  similar  discrenancies 
with  the  elementary  rigid- ideally  plastic  analysis.  Parkes 
showed  that  a  simple  correction  made  for  a  strain  rate  ef¬ 
fect  alone  led  to  good  agreement.  Both  Mentel  and  Parkes 
used  a  value  of  limit  moment  which  was  constant  in  any  one 
case,  but  multiplied  by  a  factor  to  take  account  of  rate  of 
strain.  While  their  method  of  correction  gives  the  right 
magnitudes  of  corrections  needed  in  the  cases  studied  its 
general  validity  seems  highly  uncertain. 

It  therefore  seems  desirable  to  extend  the  theoretical 
treatment  of  the  problem  to  include  strain-hardening  and 
strain-rate  effects.  The  present  report  is  an  attempt  to 
make  this  extension.  Theoretical  analyses  are  presented  for 
beams  which  exhibit  either  strain-hardening,  or  strain-rate 
sensitivity,  separately.  The  influence  of  these  effects  on 
the  permanent  plastic  deformation  is  studied.  In  addition, 
expressions  are  derived  for  the  magnitudes  of  the  plastic 
strains  arid  strain-rates  which  occur  in  the  beam,  and  numer¬ 
ical  results  are  presented  which  apply  to  the  beams  used  in 
the  experiments  [3,  4], 

2 .  General  Analysis . 

2.1  Stress-ftrain  Relations.  Two  stress-strain  relations 
for  the  plastic  range  are  considered; 

d  =  d^  +  ^  e  ,  a  >  a  (1) 

°  n  o 

e  =  D  (d  -  0o)F,  0  > 


(2) 
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Relation  (1 )  introduces  linear  strain-hardening  as  illus¬ 
trated  in  Fjg.  2.  Note  that  the  constant  n  can  be  inter¬ 
preted  as  the  ratio  of  the  slopes  of  the  elastic  and  plastic 
portions  of  the  stress-strain  curve.  Relation  (2)  introduce' 
a  strain-rate  effect  but  no  strain  hardening.  This  particula 
1  or  if;  of  the  stress-strain  relation  is  suggested  by  the  ex¬ 
perimental  data  of  Nadai  and  Manjoine  on  the  dependence  r.f 
the  yield  stress  on  the  rate  of  straining  [5,  S] .  By  suit¬ 
able  choice  of  D  and  p,  relation  (2)  can  be  made  to  agree 

• 

closely  with  experimental  a -e  curves.  This  is  shown  in  Fig. 
5,  which  compares  the  experimental  c-e  curve  for  mild  steel 
with  relation  (2). 

In  this  analysis  elastic  strains  are  neglected;  the 
material  is  assumed  rigid  when  the  stress  is  less  than  the 
static  yield  stress  cQ.  According  to  references  [l]  and  [2] 
the  criterion  for  the  validity  of  this  assumption  is  that  the 
beam's  initial  kinetic  energy  should  greatly  exceed  the  maxi¬ 
mum  elastic  strain  energy  which  can  be  stored  in  the  beam. 

The  criterion  is  satisfied  for  all  applications  which  we 
make  of  the  present  theory. 

In  order  to  develop  moment -curvature  relations  from 
the  stress-strain  relations,  the  usual  assumptions  of  beam 
theory  are  made.  Thus  we  assume  that  plane  cross-sections 
normal  to  the  centre-line  remain  plane  and  normal  to  the 
centre-line  during  bending;  that  the  cross-section  has  two 
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perpendicular  axes  of  symmetry  and  bending  takes  place  about 
one  of  them;  that  the  material  has  the  same  properties  in 
tension  and  compression,  etc.  Tnen  the  strain  and  strain 
rate  at  any  point  of  a  cross-section  are 


£  =  XZ 


£  =  XZ 


where  z  is  the  distance  from  the  neutral  axis.  With  this 
result,  and  the  stress-strain  relations,  the  bending  moment 

r  c 

y.  =  zd  ( z  .HA 
/  -c 

at  any  section  can  be  calculated  and  hence  the  moment -curvature 
relations  are  obtained.  For  the  stress-strain  relations  1) 
and  (2)  the  corresponding  moment-curvature  relations  are, 
respectively 


m  =  (l  -  \x) 


V  n  V 


X  =  j<  (M  -  y  )P,  M  >  :•!. 


wr.ere 
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o  ^  u 
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The  stress-strain  relations  (1)  and  (2 )  can  be  regarded 
as  special  cases  of  a  more  general  relation.  Malvern  [7]  has 
suggested  the  form  for  the  d-e-e  relation 

e--=g(d-f(e))  (7) 

i 

Hj 

where  d  =  f(e)  is  the  static  stress-strain  relation.  Thus, 
the  plastic  strain  rate  is  a  function  of  the  "excess  stress' 
d  -  f(e).  If  elastic  strains  are  neglected  the  term  c/E  can 
be  deleted.  The  relation 

•  „  P 

e  =  D(d  -  d^  -  i  e ) 

0  n 

is  a  special  form  of  (5))  and  includes  the  relations  (1)  and 
(2)  as  special  cases.  The  corresponding  moment- curvature 
relation  is 

h  =  D'(M  -  Mq(1  +  Xh))p 

2.2  Equations  of  Motion.  While  plastic  deformation  of  an 
ideally-plastic  beam  is  concentrated  at  certain  points 
(plastic  hinges),  the  regions  of  plastic  flow  in  a  strain¬ 
hardening  or  rate-sensitive  material  can  be  expected  to  be 
of  finite  length.  This,  plus  the  fact  that  the  position  and 
extent  of  the  plastic  regions  are  not  a  priori  known,  would 
greatly  complicate  the  analysis,  were  it  not  for  experimental 
evidence  which  suggests  that  some  simplifying  approximations 
are  permissible.  Ercm  the  tests  carried  cut  by  Mentel,  Green, 
and  Symonds  [4],  it  appeal's  that: 
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(1)  plastic  flow  occurs  only  in  a  region  adjacent  to” 
the  root  of  the  beam, 

(2)  the  length  of  the  plastic  region  is  short  compared 
with  the  length  of  the  beam, 

(p)  bending  takes  place  in  one  sense  only;  that  is, 
the  bending  moment  in  the  plastic  region  always 
has  the  same  sign. 

Thus,  'with  the  5  inch  long  beam  -specimens  tested  by 
Mentel,  Green,  and  Syr.onds,  only  that  part  cf  the  beam  within 
about  2/4  inch  of  the  root  was  permanently  bent  after  impact. 
Outside  this  region  the  beam  specimens  appeared  to  remain 
straight  and  unbent. 

The  above  experimental  evidence  may  not,  at  first  sight, 
seem  to  be  in  harmony  with  the  analyses  of  Lee  and  bymoncs 
[1],  or  Synonds  and  Leth  12],  who  show  that  two  plastic 
hinges,  one  at  the  root  ana  one  in  the  interior,  are  to  be 
expected,  and  that  the  senses  of  rotation  at  each  hinge  are 
opposite  to  each  other.  However,  Mentel  ” 3 ]  has  shown  that 
the  effect  of  the  interior  hinge  is  often  negligibly  small, 
especially  if  K  is  large.  Thus  the  rigid-ideaily-piastic 
theory  a. Iso  indicates  that  plastic  flow  is  essentially  con¬ 
fined  to  the  root  of  the  beam. 


In  view  of  tnis  evidence  ve  assume 
the  beam  after  impact  is  approximately  a 
about  the  root,  and  ve  assume  that,  for 
ing  the  equations  of  motion,  the  region 


that  the  motion  of 


rig id ’body 

rotation 

.he  purpose 

of  deriv 

of  clastic  flew  can 


oil -2b 


i r 
x  ^ 


be  ade  ;  .lately  represented  by  a  point  hinge.  Then,  treating 
the  entire  beam  as  a  rigid  bcdy,  and  referring  to  Fig.  3, 
we  find  that  the  equations  of  linear  and  angular  momentum 


are 


‘m 


=  -  (-  t/  +  G)  S  = 

2 


mb/ 

2 


(1  -  2K )  0 


( £ 


M 

m 


(i  n>p3  + 
a 


Gl)  G  =  - 


(1 


3K)  0 


(7) 


2.3  Initial  Conditions.  It  will  be  necessary  to  know  ini- 
tial  values  cf  0  and  0. 

Clearly,  the  value  cf  0  at  the  instant  of  impact  is 
0(0)  =  0.  To  find  the  initial  value  cf  Q  we  use  the  fact 
that  the  angular  momentum  about  the  root  suffers  no  impul¬ 
sive  change  at  the  instant  cf  impact.  This  is  sc  because 
the  bending  moment  Mm  is  limited  to  finite  values,  while 
the  impulsive  shear  force  Q  ,  although  possibly  infinite  at 
the  instant  of  contact,  produces  no  torque  about  the  root 
of  the  beam.  Hence  no  impulsive  torque  acts  at  the  instant 
of  impact.  Equating  angular  momenta  of  the  beam  about  its 
root  before  and  after  impact  gives 


mt2V 

2 


rJM  = 


3 


GlQ(C) 


whence 


- 

l  2i-6K 


0(0} 
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2.4  Cerent -angle  Relation  across  Plastic  hcglcn.  At  a 
particular  instant,  let  plastic  flow  be  confined  tc  the  re¬ 
gion  3  ,<:  x^  xR,  as  shewn  in  rig.  x^  is,  of  course  a 
function  of  tir.e.  For  x  >  xn,  the  beam  is  rigid.  If  is 
the  largest  value  of  x^  then  for  x  >  x0  the  bean  remains 
straight.  The  slope  at  any  point  of  the  bear,  is  given  by 

'X  rX 


O  = 


& 


do  = 


x  ax 


2a. 

5x 


In  particular, 


'*0 


hi 


xax  = 


X  CLX  "r 


0/ 


0; 


■*o 


x  dx 


x_ 


'  ^  J 
\  /  t 


Differentiation  of  (9)  gives  the  angular  velocity  of  the 
straight  part  cf  the  'Dear, 

.  ^x_  . 


9  = 


^  . 


xax  = 


XCLX 


;io) 


Oo 


since  x  =  0  for  x  >  x^. 

If  rotatory  inertia  is  neglected,  the  shear  force  and 
bending  moment  satisfy  the  equations  cf  notion 


= 

3x 


3L  = 

dx 


/  -  \ 

\  --  / 


ii’cw  deflections  near  the  root  are  snail,  and  it  is  reason¬ 
able  to  suppose  that  accelerations  near  the  root  are  small 


are 


Cll  -28 


12 


also.  Thus,  we  assume  that  y  is  so  small  within  the  plastic 
region  that  the  term  my  may  be  neglected  in  comparison  with 


Then 

Sx 


lit  =  0 

3x 


whence 


Q  -  Q 


m 


for  0  ^  x  s<  x  .  Further, 


3M  = 


dx 


6m 


h  =  -  Qmx  +  C(t) 

for  0  N<  x  ^  Xjj.  At  the  boundary  of  the  plastic  region  M=MQ ; 
that  is  M-Mo  when  x=xn,  Hence 


M  =  M.  +  (x  -  x)  Q 

II 


m 


(12) 


Putting  x=0,  we  find 


K  ~ 
m 


+  x  Q 
o  n  hn 


or 


*n 


M  -  M 

=  Urn  ;'Q 


Q 


m 


(13) 


From  (12)  we  note  that,  within  the  plastic  region, 

dx  =  -  da 

Substitution  of  this  result  into  (10)  leads  to 


(M=Ko  . 


Qm  x=C 


Q  I..=H  i 
m  1  ‘  1  hr 


ndM 


(14) 


(15) 
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If  x  is  a  function  of  M  only,  then  (15)  can  be  evaluated  to 

i 

give  a  relation  betveen  9,  if  ,  and  This  relation,  to- 

i.*  u. 

gether  with  the  two  equations  of  notion,  (6)  and  (7),  is 

sufficient  for  an  analysis  of  the  notion  of  the  bean. 

From  ( 9 )  and  ( 14) , 

- 

1  !  ^X=Xn  ,'X=X 

9  =  — '  xdM  ^  0  xdM  . 

'  x=Oj  x=xn' 


If  x  is  an  increasing  function  of  time,  then  x=0  for  x  >  x 
n  &  ’  n 

and  hence 


(16) 


This  integral  can  be  evaluated  if  x  is  a  function  of'  M  only. 
Again,  this  relation,  plus  the  two  equations  of  motion,  (6) 
and  (7),  are  sufficient  for  an  analysis  of  the  motion  of  the 
beam. 


To  investigate  whether"  ”  is  an  increasing  or  decrees- 

A  1 

ing  function  of  time,  consider 


_  -  I -6K \ 
-  -o  v - ) 


mf2(l-2K)9 


(17) 


which  is  obtained  by  substituting  for  ana  3  from  (6)  ana 

*• 

(7)  into  (13).  Note  that  G  is  a  negative  quantity. 

For  a  strain  hardening  beam,  will  increase  with  9 

•  * 

ana  with  time,  and  hence  the  deceleration  9  •.•rill  also  increase 
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in  absolute  value  with  time.  Then  it  follows  from  (17) 
that  xn  is  an  increasing  function  of  time. 

For  a  rate-sensitive  beam,  velocities  and  strain-rates 

can  be  expected  to  decrease  with  time,  and  hence  so  will  M^. 

•  • 

Therefore  the  deceleration  9  will  decrease  in  absolute  value 
with  time.  Hence,  from  (17)  ^  is  a  decreasing  function  of 
time  in  this  case. 

It  follows  that  formula  (16;  may  be  used  for  a  strain- 
hardening  beam.  For  a  beam  which  exhibits  both  strain- 
hardening  and  strain  rate  effects,  however,  it  is  not  clear 
how  xn  changes  with  time,  and  hence  it  is  uncertain  whether 
formula  (16)  applies. 


3*  Analysis  for  Strain-Hardening  Beam. 

3.1  Calculation  of  Damage  Angle  Of.  Here  the  moment-curva¬ 
ture  relation  is  (3),  namely 

M  =  MQ(i  •-  \x ) ,  M  >  Mq 


or 


x 


M  -  M 


a 


\M0 

Substituting  this  into  (16)  and  integrating,  we  obtain  the 
relation 


9 


<«-  -  Kn>2 

2X  V‘o 


(18) 


Elimination  of  y  between  (18)  and  the  equation  of  motion  (6) 
and  (7)  leads  to  the  equations 


Cll-23 


15 


G  =  (£l6Kj  i_  (3-lj: 

3+6K  2X  S 


3 


(19) 


e  =  - 


6M, 


ml3(2+6K) 


where  3  =  M  /K  »  Solving  the  first  of  equations (19 )  for  S 

v 


we  obtain 


S  =  (1  :  uG>  !  V)'  (1  +  aGT  -  1 


wnere 


V 


a  _  X  (3+6K}  _  El  ,3+6K, 
2+6K  nM0^  2+6K 


(20) 


Substitution  of  (20)  into  the  second  of  equations  (19)  gives 

r  , _  i 

••  Am  1  ,i 

9  =  - 


mE3(2+6K) 


12— _  '(1  +  a0;  I  J(  1  +  aGT  -  1  (  (21) 


ns  G  increases  so  does  M^,  because  of  the  strain  hardening 
and  consequently  the  angular  deceleration  of  the  beam  in¬ 
creases  also.  This  shows  that  the  positive  sign  must  be 
chosen  in  (21), 

Now 

••  jA  •  iA  j  *9 

(22) 


29  =  2  =  2Q  =  4-  (q2) 

at  d9  dG 


and  hence,  from  (21)  and  (22), 

r 

12M 


d(02)  = 


xl 


9 


^1  +  aQ  +  ;(1  +  a9r  -  1  ,vdG 


2 


ml- ( 2+6K ) 


\  I 


(23) 
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The  damage  angle  ©^  is  found  by  integrating  (23).  The 
limits  of  integration  are:  at  the  beginning  of  motion  9=0  and 

r\  V  /  3  -t  ^  Y  “ 

y  ~  r  (oT^7);  the  motion  ceases  when  9  =0,  at  which  time  9  =  C*. . 
Hence 


Q=0 


d(a2)  = 


-12M, 


9, 


‘O  ;  I  •  /  o 

- - —  ,;1  +  cr©  +  w(l  +  a©r  -  1  ' dQ 

m^(2+6K)  •  1 


or 


n  +  6K)3  = 


m£IV2 

12 nM2  (2  +  6K)2 
o 


0 


aQ^ 


<il  +  K  +  \/(l  +  T;2  -  1  'dS 


0 

1) 

r 

i 

0 


\  1  +  z  +  i/(l  +  K)  -  1  \&l  (24) 


where 


^  =  ag 


n  =  aGp  =  (3 — +„6li )  q 

1  f  il  2  -  6K'  yf 


Now 


n 

'0 


1  +  K  +  Vd  +  O'  -  l  ;d£ 

-  Js  Jon  _1_  „2 
2 


(25) 


=  -  FCn) 

2 


<! 2 rj  +  if  +  (1+n)  \l(l+r)  )d**l  -  ^n(l+T)+  ./(l+n)^-l) 

(26) 


Then  (24)  becomes 


?<n)  =  afcjpd  a  *  6k)3 

6P.K2  (2  -  6  y.)2 

0 

=  £  u  -  2!'-)3 
8  2  2 
nM  (1  +  3K) 
o  ^ 


and  (25)  nay  be  vritten, 


using  (27) 

a  ~  n 

-f  an 

r*  “  ?TnT 


where 


a 

-  n 


3  ^y2  (i  +  2  k)2 
S’  K0  (1  -“??0  ' 


'27^ 


(28) 

(29) 


is  the  damage  angle  which  would  be  obtained  if  the  beam  were 
Ideally  elastic.  Relation  (29)  is  derived  In  reference  [4], 
Equations  (27)  and  (28)  give  the  relation  between 


9 ~/Q  and  the  non-dimensional  velocity  S.EI7  (1  +  2h)E 

nir  (1  +  “*Kj 

o 

relation  is  plotted  In  figure  6,  together  with  experimentally 
determined  points  taken  from  the  data  of  reference  [3]* 

In  order  to  plot  the  experimental  points  an  estimate  ci 
the  value  of  n  is  needed,  which  is  obtained  as  follows.  V.’e  take 
the  slope  of  the  post-yield  portion  of  the  stress-strain  curve 
to  be  approximately 


6 


TJL 


~  o 


0 


E  ULT 


where  d  m  is  the  ultimate  stress  of  the  metal  and  e  the 
ULT  ULT 

corresponding  fractional  elongation.  Then 


n 


n  e 


ULT 


-  d 


o 


C 1 1  -  2  8 


lb 


For  tne  mild  steel  beams  used  in  the  experiments,  d  1  5 6,000 

U  iji 

psi.,  <3  21  33*000  psi.,  e  1  0.1,  whence  n  1  130.  This  value 

was  used  in  plotting  the  experimental  points  of  figure  6. 

The  non-dimensional  velocity  ffiSI’/f.  XX  1  can  ba 

nM2  (1  +  3K)2 
o 

interpreted  in  terms  of  the  ratio  R  of  the  beams  initial  kinetic 
energy  to  the  maximum  elastic  strain  energy  v/hich  can  be  stored 
within  the  beam.  (Recall  that  it  is  this  ratio  which  determines 
whether  elastic  strains  may  be  neglected.)  The  expression 


for  R  is 


\  ml’2  +  \  GV2 

R  = - “FI - 

1  o 

2  El 


=  ( i  +  k) 

% 


and  hence 


a  +  2*)' 


i 


mETV2  <i  "  2KT  _  B  _  j _ 

O  +  3K)2  n  '  l(1  +K)(1  +3k)2J 


The  quantity 


!  (1  +  2K)3 

2 

_(1  +  K )  (1  +  3K) 


is  plotted  in  figure  7*  and  it 


can  be  seen  that  this  quantity  is  almost  constant  over  most  of 
the  practical  range  of  K.  Thus  the  reduction  in  damage  angle 
due  to  strain-hardening  is  chiefly  a  function  of  R/n,  the  para¬ 
meter  K  having  only  a  small  influence. 
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3,2,  faxculaticn  of  .-iaxlmmn  Plastic  Strain 

The  maximum  strain  occurs  at  the  root  of  ti¬ 
the  rending  moment  is  greatest. 

From  (3)  i  v; 

x  =  ^  ( rr~  ~  1 ) 

A  *'*0 


^  »*«  *  *0  0 


and  hence  the  curvature  at  the  root  of  the  beam  is 

h  =  f  !rf  -  1)  =  — 

m  X  ..0  X 


rr  ^ 

—  Qg 


u/  ( 1  •+■  a  3  )  ^ 


\ 


using  (20),  The  strain  at  the  extreme  fibre,  at  the  it.  o'- 


then 


x  C  =  £  .•  a9  +  J  (1  +  a9 ) 
m  X 


> 


and  the  final  strain  at  the  end  of  cne  deformation  is 


t  'a9f 


I  2 

7(1  +  a9j 


-  1 


M  o 


n 


j  (If 


Equations  (27)  and  (31)  give  the  relation  between  the 
maximum  strain  and  the  non-dimensional  impact  velocity.  This 
relation  is  plotted  in  figure  8  using  values  of  the  parameters 
aupropriate  to  the  mild-steel  beams  which  were  used  in  the 
experiments  reported  in  reference  [3]- 


c*j 
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3.3. 


Calculation  of  Orain-Ra  tes  : 
From  (30)  and  (31), 
nNLc 


c  = 
m 


o 

El 


[a©  +  \/(l  +  a9  f  -  1  1 


and  differentiation  of  this  result  with  respect  to  time  gives 


nM  c 

4»--ir  [1  + 


]  a0 


(32' 


•/(l  +  a©)2  -  1 

Note  that  as  9  tends  to  zero  and  0  tends  to  the  initial  (finite) 

angular  velocity,  then,  according  to  (32),  e  tends  to  infinity, 

m 

Thus  the  present  theory  predicts  that  the  initial  strain-rate 
is  infinite. 

The  moan  strain-rate,  i.e.  the  time  average  strain-rat;: 
at  the  root  of  the  beam  is 

f. 

edt 


e 


m 


MEAN 


(33) 


dt  j  dt 


where  the  integration  extends  over  the  time  during  wh ich  plastic 

flow  occurs.  To  find  I dt,  return  to  equation  (23)  and  integrate 

•  •  • 
it  between  the  limits  9  =  9,  9  =  9,  and  the  final  values  9=0, 


9  =  0^.  Thus 


9=0 

I  *2 

d(e  )  = 

/ 

9=9  m 

whence 


-12M. 


s=ef 


lh2  -  6K)  4 


j  Vjl  +  cc9  +  /(I  +  a©)  -  1  fdO 


•2 

G  = 


6nNu 


o 


mEl£(3  +  6X) 


[F  ( r) )  -  F(a© )  ] 


using  (26)  and  (25) 


/if ter  solving  for  9  =  d9/dt,  we  can  obtain 


"11-2P 


It  = 


cnM2 


C F( n  )  -  F(a0 ) ] * 


The  length  of  tine  luring  which  plastic  flow  occurs  then  is 


dt  = 


bnM2 

o 


,  9=e^ 

T.i  i  ;  -4 

L  1  [F(r)  -  F(a© ) ]  d© 


mntM  (2  +  6K)2  ^ 
6FT  (3  +  6K) 


[f ( n )  -  F(Ol 


(2*v5 


/  K=o 

on  substituting  £  =  a©  from  (25). 

The  integral  In  (3*+)  has  to  be  evaluated  numerically. 
In  the  above  form,  however,  the  integral  is  not  suitable  for 
numerical  treatment  because  of  the  singularity  at  the  upper 
limit.  This  singularity  can  be  removed  by  integrating  once  by 
parts  ana  then  making  the  substitutions 


£=!+?-  y(l  +  u2  -  1  ( 


p  =  1  ♦  1  -  y  Cl  +  n)2  -  1 


This  maneuver  reduces  (3*+)  to 


(35) 


m^2V(3  +  6K)  J 


dt  =  1  - 


[1  -  H(^  V  dC 
L  hTpTj 


(36) 


while  the  substitutions  (35)  change  (26)  and  (27)  to 


H(P  >  =1 


IV2  (1  +  2K)3 


n!£  (1  +  3k)‘ 


(37) 


cu-aj 


22 


where 

H(p)  z  -»  -trip2  -  1  (38) 

P 

finally  the  time-average  strain-rate  is  found  using  (31  ),  (33 ), 
(36),  (37),  and  (36),  as 


n f  i  \ 


1  pti  _ _ 


(39) 


MIAN  l  ajE!  (l+3K)2  '  f  2  p  2  2 ■  *  h  JJTcTXT” 

-*  ,i+p  ^p  -p  ,  ■;!-  r[i-  ^(“jJ  uC,> 

Equations  (37)  and  (39)  give  the  relation  betv/oen  non- 
dimensional  impact  velocity  and  mean  strain-rate.  This  is 
plotted  in  figure  9,  again  using  values  of  the  parameters  appro¬ 
priate  to  the  mild-steel  beams  which  were  used  in  the  experiments 
reported  in  reference  [3]*  We  note  that  the  mean  strain-rate 
depends  significantly  on  K. 


3 . 4- .  Length  of  Plastic  Redon 

From  the  equations  of  motion  (6)  and  (7)  we  obtain 


0  = 

1  ^2  +  6%J 


and  substitution  of  this  result  into  (13 )  leads  to 


=  (2  1  6K)  £. . 1 
t  o  +  or  s 

where  S  =  Mm/MQ.  Substitution  for  S  from  (20)  gives 

-f  -  “9  +  l/(l  +  a@)2  -  1) 

The  maximum  value  of  xn  is  given  by 


(40) 
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ir.axffx!)  =  (f — 6K > ( “  a8f  -  /(l  +  cSf)2  -  1 

=  (|  ^  |%)(-  n  +  .'a  *  n)2  -  i)  (4i) 

X 

Vhe  value  of  r.ax(_|i)  is  plotted  against  the  non-dinenslonal 

impact  velocity  in  figure  10. 

Since  it  is  a  basic  assumption  of  the  present  analysis 

that  x  /t  is  small,  figure  10  shows  that  there  is  an  upper  limit 
n 

on  the  impact  velocities  for  which  the  present  theory  is  valid* 

The  requirement  x  /t  «  1  will  be  satisfied  provided,  say, 

n 

smi  iLizai  <  .05 

nK2  (1+3  K)2 

Recalling  from  the  liscussion  of  section  3*1  that 

raEIV2  (1  +  2K)3  ~  ^  R 
nM^  (1  +  3K)2  n 

then  the  requirement  can  be  stated 

R  <  .06n 

For  n  =  130,  this  becomes  R  <  8,  roughly.  On  the  other  hand, 
there  is  the  conflicting  requirement  R  >>  1  which  must  be  satis¬ 
fied  in  order  that  elastic  strains  may  legitimately  be  neglected. 
Therefore  it  appears  that  the  present  analysis  is  valid  only 
for  a  limited  range  of  R,  in  the  neighborhood  of  R  =  8. 


a  •  nr.alvsls  for  Rate-lens J  tlve  Beam 

■+.1 .  J  -Iculntion  of  damage  angle  Q* 

Hero  we  use  the  mcn.^nt-curvature  relation  (4),  which 
introduces  a  strain-rate  efi'ect, 


x  =  D»  (M  -  M0)‘ 


=  4  CM-  -  DP 

O  » 

'  '0 

where,  for  a  beam  cf  rectangular  section, 


(4 


P  =  cD'MP  =  Ddp  (£P.-..t.-l)P 
o  o  2p 

Vfc  note  that,  in  te^ms  of  p5  the  stress-strain  relation  (2)  can 

be  written  .  ~  n 

e  =  |i(— (-2.  -  1)P 

From  (42)  and  (15)  we  find  that  the  moment-angle  relation 
across  the  clastic  zone  is 


M 


m 


c  -  JL 


©  =  rf-  1  x  (M)dM 


J 


M. 


PM, 


*V 


o  ( ' m  n  xP 
Qm(p  +  l)c  V 


(43 ) 


"o 

Elimination  of  Qm  between  the  equations  of  motion  (6),  (7)  and 
(43)  leads  to 


9  =  gj  .  IS  -  liH 

(p  r  l)c  j  +  6K  S 


-  bM 

9  =  — - 2- 


m  &  ( 2  +  6  ^ 


( 4if ) 


(45) 


where  S  =  M  /M  .  From  (44), 

-mo  ? 
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Substitution  from  (44),  (45),  and  (46)  into  (47)  yields 

d©  -  -  -  m3  1°  (2  +  6K) 

6(p  +  1)c2Mq  (3  +  6K)2 

(48) 

The  damage  angle  ©^  is  obtained  by  integrating  (48).  The  limits 
of  integration  are  obtained  from  (44).  At  the  beginning  of  the 
motion  9=0  and  9  =  ^  ^2+6K^ •  Then,  using  (44)  we  find  that  SQ, 
the  initial  value  of  3,  is  given  by 


(S0  -  1)  2 

TT-,  ^  =  U  (^SS) 


(49) 


(p  +  1)S0  pfc  2  +  6K  '  " 

•  • 

The  motion  ceases  when  9=0.  Alien  9  =  0,  then  3=1,  from  (44), 
Integration  of  (48)  between  the  limits  9  =  0,  S  =  S0  and  ©  =  0f, 


S  =  1  gives 

n  _  -  (2+ 6K)3 

y  '  ' r” 

6(p  +  l)c2Mn  (3+  6K)2 


.3  -  1 


2p+l 


1  (3-1 


2p^2 

'  ■  ■»  / 


=  0.  _ £ _  fiib  tZ.+  .  ! 

°  (p  +  i)  yr  A  +  6 kj  ■  r : 


2p+l 


3+6 


(P+D  j 

-.1.  ISz I 

(p+i)  ; 


>  ( 

s  i 

2p+2 ; 


(50) 


where 


i— i!  oj 
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is  the  damage  tr.gle  which  wculd  be  obtained  if  the  beam  vv  re 
ideally  p]astio.  -hjosti tution  of  (49)  into  (50)  gives 


«  f 

/  =  2(P-1) 


(s0  -  1) 


2  p+2 


u  2p+l  2p+2 

.  ,  (S  -  1)  H _ 1  (S-l)  p 

.  |  d3  TFTT  jr  j 

l  1  ( 51 ) 


1  dS 


By  integrating  the  second  term  in  the  integral  by  parts  ve 
obtain  ~ 


2(P  +  ds; 


90  "  3S 


0  i(s0  -  i)2p*2  Ji 


0  2p+l 

,  IQ  ■  -  U 


dS 


(  52) 


If  2p+l  is  integral,  then  the  integration  of  (52)  is 
elementary.  However  the  results  are  not  suitable  for  computation 
and  a  better  procedure  is  the  following.  By  introducing  a  new 


dummy  variable 


(52)  is  reduced  to 


S  -  S 


_  o 


So  -  1 


,1 


7T  =  tI“  +  '2~^t c  i  (1  -  u)"1'  [1  -  u(— - - )]  du 


3SC 


3  S0 


2p+l 


- 


-3 


Expansion  of  [l  -  u( 


0 

So  -  1\  -r3 


)  ]  by  the  binomial  theorem  then  gives 


if  =  _2_  + 

e0  3S0  is 


(l-u)2p+1  2  (ni-l)(j°_i)n  un(iu  (53 

n=0  2  S, 


■s0-l.n 
3o 


For  0  <_ u  <_1  ,  the  series  above  converges  uniformly  for 

<  SQ  < co  .  This  range  includes  all  physically  possible  values 
of  S0.  Interchange  of  the  order  of  summation  and  integration  in 
(53)  yields 


■u-id 


9  3s 
0  J  0 


O  n  00 

-4 ILJL-i  7 

3so  n=C 


!  (1-u) 


P  +  1 


*,  *  ^  1  ‘ 

L*.  wi 


Using  the  known  result 


'0 


u^Ui  -  uT-Eu  =  vvrou 

r( a  +  b) 


we  obtain 


2 

3S, 


®  -Ln^^.u  (fcJ:}n  n.uiij.r£2p.iii 

^  o  n=C  ^  i Up  +  n  ^  j ) 


2p  +  2 


(  oo 

i  J  i  t  i  2  r(nt3jr(?.p~u3  2  ( 
°0  '•  6  n=l  r(2p  -  n  +  3  )  ^ 


-  1  n'< 
-V— )  > 


(5^) 


The  series  in  (^h)  converges  rapidly  and  is  suitable 
for  numerical  calculation.  Note  p  need  not  be  integral.  Equa¬ 
tions  (4-9)  and  ( 5M  are  parametric  equations  relating  the  damage 
angle  9^  to  the  impact  velocity  V,  through  the  parameter  So<> 
This  relationship  is  shown  in  figure  11.  Experimental  results 
from  reference  [3]  are  also  shown  for  comparison.  In  plotting 
the  experimental  points  the  same  value  of  (J  was  used  as  was 
used  in  drav/ing  figure  5 • 

'tie  note  that  the  theory  predicts  a  sharp  drop  in  the 
damage  angle  due  to  the  strain-rate  effect.  Over  most  of  the 
range  of  velocities  the  damage  angle  is  about  half  that  pre- 
dieted  by  the  rigid-iaeally  plastic  theory,  which  is  in  rough- 
agreement  with  experimental  results.  The  experimental  points 
all  lie  below  the  thcorerical  curves,  and  this  might  be  attrib¬ 
uted  to  the  effect  of  strain-hardening  which  we  have  not 


cii->e 


^8 


considered  here. 

Three  parameters  appear  in  the  stress-strain  relation, 

namely  p,  K0  (or  d0)i  and  j3  .  In  the  curve  of  damage  angle  vs. 

non-dimensional  impact  velocity  the  parameter  MQ  enters  only 

in  the  scale  factor  of  the  ordinate  S^/0O,  while  the  parameter 

3  enters  only  in  the  scale  factor  of  the  abscissa  YSL- 

a  12  2+6K 

Only  the  exponent  p  affects  the  shape  of  the  curve.  Figure  11 
Indicates  that  the  curve  is  rather  insensitive  to  changes  in  th 
value  of  p. 

Changes  in  the  value  of  {3  would  have  the  effect  of 
shifting  plotted  experimental  points  to  the  right  or  left  in 
figure  11.  Ginee  the  curves  are  rather  flat  over  most  of  the 
range,  changes  in  the  value  of  (3  would  do  little  to  move  experi 
mental  points  closer  to  or  farther  from  the  theoretical  curves. 
Thus  the  agreement  between  theory  and  experiment  v.’culd  be 
influenced  little  by  inaccuracies  in  the  value  of  j3  ,  The  ordi¬ 
nate  90/©.p,  however  involves  MQ(or  dQ)  directly.  Changes  in 
M0  would  have  the  effect  of  moving  plotted  experimental  points 
upward  or  downward,  and  this  would  strongly  affect  any  agree¬ 
ment  or  disagreement  between  theory  and  experiment. 

Therefore,  in  the  present  theory,  the  most  critical 
parameter  is  Mq  (or  the  others  seem  to  be  rather 

unimportant . 


I»-cO 


~Q 
—  / 


*♦  .3  .  :alculv.icn  of  Maximum  Plastic  Strain 
The  strain  is  calculated  from 


e  =  t  edt 


'  crc  \ 

7  J  / 


.-it  tie  roct  of  the  beam 


E  =  X  C  =  -  1)P 


m  □ 


from  (42 ) ,  This  is  the  greatest  strain-rate  within  the  tlastic 

O 

regicn.  Then  we  have 

t 

M  p 

=  ?■;  (J  - 1)  at 

'  o 
3 


=  ?\  (S-l)P^ 


0 


A 


since  at  =  dG/S<  Hence,  using  (44 )  and  (48), 

.3^ 


,4.2 

-  o  u 


rr *  A  V 

^  ^"0 


o  (3  +  6X; 


9  2d 

Ifbt  \  I  (£  -  1)  1 


P  ^  1) 


2dt1  ! 


dS 


(56) 


The  procedure  used  to  evaluate  (56)  is  essentially  the  same  as 
the  procedure  used  to  evaluate  ( 51 ) •  Going  directly  to  the  final 
result, 


=  BcVt  H 


1  P  >/ 

‘  *0 


A-.,3  ' 

JS4.L,  42- 

/  / 


U 


(2  -r  6K 


2!  O 
v° 


-  1) 


1 


o 


CD 


n=0 


nm mi ! 

“  -  2) 


(57) 


(57) has  been  evaluated  using  values  of  the  parameters  appropriate 
t''  the  mild-steel  beams  used  in  the  experiments  of  Mentel< 
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Green ,  and  Synrr.ds  [3],  [4-],  Also,  the  values  of  p  and  [3  used 
in  plotting  figure  5  were  adopted.  The  results  are  snc.’Ti  in 
figure  12. 


Note  that  the  strain  calculated  here  is  the  strain  at 
the  root  of  the  beam,  which  is  the  maximum  strain.  The  average 
strain  in  the  plastically  deformed  section  of  the  beam  is,  of 
course,  less  than  this. 

It  is  seen  from  figure  12  that  the  curves  of  e  vs.  v 

b  m 

are  very  closely  straight  lines  on  a  logarithmic  scale.  This 

suggests  that  the  e  -  V  relation  can  be  acproximated  by  a 

m 

power  !  w.  Using  only  the  first  term  of  the  series  (57)  we  have 


-  mc.V  jQ 
m  -  24M 


+  6K)3 


Up  +  3 


0  ,  2  (3  -  1. 

0  (2  +  6K f  0 


( 53 ) 


From  (49) 


<S0  -  1)  =  [-fr-i-ps  (Hrl)2]P+1  s?*; 


$1‘ 


(59) 


Use  the  value  p  =  5»  Then  substitution  of  (59)  into  (58)  gives 


m  — 


m«h/y/6,l/3  n  +  ac  )6/3  v11/6 

(2  +  6K)5/3  S  1/6 
0 


2.49Mr 


(6C) 


The  range  of  interest  of  SQ  is  1.6  <  SQ  £  2.0  roughly.  In  this 
range  S  varies  by  only  b%.  Hence,  putting  3^  ^  ^ 

=  1.103  in  (60)  we  obtain 


m  — 


.S3B.l67|,33n  „  6KF-67  1.83 


2.74K0(2  +  6K)1’67 


(61) 


as  an  apcroximate  formula  for  the  maximum  strain. 


^  ^  -  -  -  V-/ 


4.3.  3 3 i dilution  of  3tr'-ln  Rates 

;-.t  the  root  of  the  beam 

•  •  M  p  p 

e  =  xc  =  ?  (-a  -  i)  =  ,6(s  -  1) 

“0 

from  (43).  The  maximum  strain-rate  occurs  at  the  beginning  of 
the  sefcrmaticn  (when  5  is  largest)  and  hence  is 


£ 


H«X 


?(30  -  1, 


■  P 


(62) 


r.quat; ons  (49)  and  (62)  give  a  parametric  relation 
between  maximum  strain-rate  .and  impact  velocity. 

«gain  an  aprroximate  formula  can  be  developed.  Substii 


tien  of  (59)  into  (62)  yields 

1 

•  p-t-l  (p  »■  1 )  V  cS 

e.-/=p-  [ - - - -  t 


;■ 


O 


0  ^  -  6K  2 


2  -  6K 


)  ] 


Taking  0=5,  and  using  for  3  its  mean  value  1.6,  ve  obtain 


.167  .83  1.67  5- 

e  2«.5P  _ S _  r3  +  6h> 

^u»X  -  1.67  ^ 


2  +  6K ' 


V 


f  '  -j  \ 

vod  j 


0 


The  mean  strain  rate,  i.e.  the  time  average  strain  rate, 
at  the  root  of  the  beam  is 


6  MLrtN 


J  m  e 


m 


(64) 


j  dt  j  dt 

where  the  integration  is  taken  over  the  interval  during  which 
plastic  deformation  occurs.  has  already  be^n  evaluated.  Tow 


dt  =  ££ 


=  1T10L4I  \  J  13-1) '  _  1  (o-l )?  1 

1  ,  62  tfh  ,3 

1  ^ 


_ _ ILL 

6 MQc  (3'6KT 


dS 


■  a  c  •- 

\  ^  y 
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by  using  (44)  and  (48).  The  integral  (65)  can  be  evaluated  by 
the  same  procedure  used  in  evaluating  (5D.  The  result  is 


m 


00 

1  +  z  CXn 
n=o 


^ o” ^  n 


n+2jnp.j 2 1 

n  +  p  +  2)  K  3o  ' 


(66) 


Substitution  of  (66)  and  (57)  into  (64)  yields 


MEAN 


„  p  (p+l)Sn  I 

vc,i±£J^2  * _ £  _! 

^2^2+6K  <S0  -  1) 


]  1  + 


x  iwm*  ¥)0 ; 


4  +  y  r(n+2)r(wi  (!o±)n 

(  a=C  ^(n  +  p  +  2)  S 


£  MAX 


1+2 


'n+2)r(2p+l 


n=0  r(2p  +  n  +  2 


\i  +  z  rintfjr-Ca-t&l 
[  n=0  rTnT  p  +  2) 


(67) 


in  view  of  (62)  and  (49). 

•  • 

In  figure  13,  numerical  results  forew>„,  and  eWT,.,.T  are 

MAX  MEAN 

shown.  We  note  that,  within  the  range  of  velocities  considered, 
£MEaN  is  very  nearly  one-half  of  e  NlAX« 


4.4,  Length  of  Plastic  Zone 
From  (40) , 


t-  “  X  +  6K;  S 


The  maximum  value  of  x  /-l  occurs  at  the  beginning  of  the  motion, 

n 

and  is  v  S  -  1 

maxCSU  -  (2_t||)  (63) 
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The  mean  value  of  x^/t,  i.e.  the  time  average  of  is  given 

by  f  x_ 


XL  I 

mean  ( -It)  =  . — 

7.  i 


(69) 


The  expression  (69)  can  be  evaluated  in  the  same  way  as  the 
expression  (64)  for  the  mean  strain-rate  was  evaluated.  Going 
directly  to  the  final  result, 


mean(Jl)  =  j,  J 

l  3+6K  Sq  f 


J  l  +  s  Clh:.3)r(l,t2J. 

!  ^  n=0  P(n  +  p  +  3) 


(£o±}  s 
+  3)  S0  * 


I  CO  c  _T  n 

,1+2  n(n+2^r(p+2)  . 

n=0  T(n  +  p  +  2)  SQ 

1  +  X2bui  f  g ia+uctegi  (!cl)n 


=  ;.-,x(-jQ)  •  | 


,-n  P(n  +  p  +  3 


f  00 

<  1+2 


2  C(n+2)r(p+2 )  (s0”1)n  , 

n=0  T(n  +  p  +  2) 

r  0  2 


(70) 


In  figure  14,  max(--Jl)  and  mean(lU)  are  plotted  against 

l 

the  non-dimensional  impact  velocity.  It  is  seen  that,  according 

to  the  present  theory,  the  mean  value  of  x  /£  is  about  1/4  over 

n 

most  of  the  range  of  velocities,  and  that  the  peak  value  of 
x  /'t  can  be  as  high  as  1/2,  Therefore  it  appears  that  the  basic 
assumption  x  /£  «  1  is  satisfied  only  in  a  very  crude  manner, 
and  consequently  the  accuracy  of  the  present  theory  is 


questionable . 
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5 .  Discussion  of  nesults 

The  foregoing  analyses  were  based  on  the  assumption  ? 

among  others,  that  the  length  of  the  plastic  zone  is  short 

compared  with  the  length  of  the  team.  All  test  results  with 

both  mild  steel  and  aluminum  alloy  specimens,  described  in 

references  [3]  and  [4],  indicate  this  assumption  to  be  correct. 

The  region  of  large  plastic  deformation  occurs  in  a  length 

approximately  C.5  in.  or  less,  so  that  experimentally  x  /Z  <_  0,1. 

However  the  analysis  developed  in  this  report,  which  is  based 

on  this  assumption,  predicts  that  under  the  conditions  realized 

in  many  of  the  tests  the  value  of  x  /Z  is  actually  so  large  as 

n 

to  invalidate  the  theory.  The  reasons  for  this  are  not  yet 
understood,  and  the  analysis  is  unsatisfactory  to  this  extent. 

The  discrepancy  may  perhaps  be  due  to  the  use  of  stress-strain 
or  stress-strain  rate  relations  which  are  unsatisfactory,  since 
it  is  by  no  means  certain  that  stress-strain  data  obtained  by 
Man jo in e  [6]  for  high  strain  rates  can  be  applied  to  problems 
of  bending  in  the  manner  done  here.  Further  development  of  the 
theory  with  a  view  to  removing  the  assumption  xn  <<  Z  is 
desirable . 

Both  the  analyses  accounting  for  strain-hardening  and 
strain-rate  sensitivity  separately,  predict  that  the  damage 
angle  is  considerably  less  than  the  angle  calculated  on  the 
basis  of  a  rigid-ideally  plastic  material.  Of  the  two  effects, 
strain-rate  sensitivity  appears  to  be  the  more  important,  and 
this  conclusion  is  strengthened  by  the  findings  of  Manjoine  [6] 
that  the  stress-strain  curve  becomes  flatter  at  higher  strain- 
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rates.  The  reduction  in  damage  angle  due  to  either  strain¬ 
hardening  or  strain-rate  sensitivity,  by  itself,  is  not  sufficient 
to  bring  the  thoeretical  analysis  into  agreement  with  experimental 
results.  From  the  magnitudes  of  the  reductions  due  to  the  two 
separate  effects,  however,  it  seems  likely  that  agreement  could 
be  obtained  by  considering  the  combination  of  the  two  effects. 

This  has  previously  been  stated  by  Mentel  [3].  Thus  it  would 
be  desirable  to  extend  the  present  analysis  so  that  strain¬ 
hardening  and  strain-rate  sensitivity  are  accounted  for  simulta¬ 
neously. 

That  an  extension  of  the  analysis  in  this  direction  is 
desirable  is  evident  also  from  the  curves  of  strain  and  strain- 
rate  vs.  impact  velocity.  Thus,  the  predicted  mean  strain-rates 
in  a  strain-hardening  beam  are  so  large  that  the  neglect  of 
strain  rate  effects  is  not  justified.  On  the  other  hand, 
predicted  strains  in  rate-sensitive  beams  are  so  large  that 
some  account  should  be  taken  of  strain-hardening,  even  though 
its  importance  diminishes  at  high  rates  of  strain. 
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FIG  8  PREDICTED  STRAINS  IN  STRAIN  -  HARDENING  MILD  STEEL 


BEAMS 
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FIG  12  PREDICTED  STRAINS  IN  EXPERIMENTAL  MILD 
STEEL  BEAMS,  ACCOUNTING  FOR 
STRAIN- RATE  SENSITIVITY 


DATA:  • 

68.0  sec.” 1 
i  -  5  in. 

c..  =  i/16  in  . . 


Cll  -  28 


FIG  13  PREDICTED  STRAIN-RATES  IN  EXPERIMENTAL 
MILD  STEEL  BEAMS,  ACCOUNTING  FOR 
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